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A Reduction Method for Higher Order Variational 
Equations of Hamiltonian Systems 

A. Aparicio Monforte and J.-A. Weil 

Abstract. Let k be a differential field and let [A] : Y' = AY be a linear 
differential system where A G Mat(n , k). We say that A is in a reduced form 
if A £ g(k) where g is the Lie algebra of [A] and k denotes the algebraic 
closure of k. We owe the existence of such reduced forms to a result due to 
Kolchin and Kovacic Kov71 . This paper is devoted to the study of reduced 
forms, of (higher order) variational equations along a particular solution of 
a complex analytical hamiltonian system X. Using a previous result |AW| . 
we will assume that the first order variational equation has an abelian Lie 
algebra so that, at first order, there are no Galoisian obstructions to Liouville 
integrability. We give a strategy to (partially) reduce the variational equations 
at order m + f if the variational equations at order m are already in a reduced 
form and their Lie algebra is abelian. Our procedure stops when we meet 
obstructions to the meromorphic integrability of X. We make strong use both 
of the lower block triangular structure of the variational equations and of the 
notion of associated Lie algebra of a linear differential system (based on the 
works of Wei and Norman in AVN63J. Obstructions to integrability appear 
when at some step we obtain a non-trivial commutator between a diagonal 
element and a nilpotent (subdiagonal) element of the associated Lie algebra. 
We use our method coupled with a reasoning on polylogarithms to give a 
new and systematic proof of the non-intcgrability of the Henon-Heiles system. 
We conjecture that our method is not only a partial reduction procedure but a 
complete reduction algorithm. In the context of complex Hamiltonian systems, 
this would mean that our method would be an effective version of the Moralcs- 
Ramis-Simo theorem. 



1. Introduction 

Let (k , ' ) be a differential field and let [A] : Y' = AY be a linear differential 
system with A G .M n (k). We say that the system is in reduced form if its matrix 
can be decomposed as A = $D j=1 onAi where on 6 k and Ai G Lie(Y' = AY), the 
Lie algebra of the differential Galois group of [A] . 
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This notion of reduced form was introduced in |Kov7l| and subsequently used 
(for instance jMS96a and MS96b ) to study the inverse problem. It has been re- 
vived, with a constructive emphasis, in [AWj . It is a powerful tool in various aspects 
of linear differential systems. The main contribution of this work lies in the context 
of Hamiltonian mechanics and Ziglin-Morales-Ramis theory MRRS07 : reduced 
forms provide a new and powerful effective method to obtain (non-)abelianity and 
integrability obstructions from higher variational differential equations. 

This article is structured in the following way. First we lay down the back- 
ground on Hamiltonian systems, differential Galois theory, integrability and Morales- 
Ramis-Simo theorem. In section [3l we define precisely the notions of reduced form 
and Wei-Norman decomposition and the link between them. Section |4] contains the 
theoretical core of this work: we focus on the application of reduced forms to the 
study of the meromorphical integrability of Hamiltonian systems. We introduce a 
reduction method for block lower triangular linear differential systems and apply it 
to higher variational equations, in particular when the Lie algebra of the diagonal 
blocks is abelian and of dimension 1. In section [5l we demonstrate the use of this 
method, coupled with our reduction algorithm for matrices in sp(2,k) [AWj by 
giving a new, effective and self-contained Galoisian non-integrability proof of the 
degenerate Henon-Heiles system f |MR99] . |MRRS07] . |MS09p which has long 
served as a key example in this field. 



2. Background 

2.1. Hamiltonian Systems. Let (M , u>) be a complex analytic symplectic 
manifold of complex dimension 2n with n G N. Since M is locally isomorphic 
to an open domain U C C 2 ", Darboux's theorem allows us to choose a set of 
local coordinates (q , p) = (qi ...q n ,pi.. .p n ) in which the symplectic form u is 

expressed as J :— 7 q . In these coordinates, given a function H £ C 2 (U) : 

U — > C (the Hamiltonian) we define a Hamiltonian system over U £ C 2n , as the 
differential equation given by the vector field Xh '■= JViJ: 

(2.1) <?i = ff(<7,p) , ft = -§§(<?, P) fari=l...n 

The Hamiltonian H is constant over the integral curves of (|2.1[) because Xh • 
H := (VH , Xh) = (VH , JVH) = 0. Therefore, integral curves lie on the energy 
levels of H. A function F : U — > C meromorphic over U is called a meromorphic 
first integral of (12.11) if it is constant over the integral curves of (|2.1I) (equivalently 
Xh ■ F — 0). Observe that the Hamiltonian is a first integral of (|2.ip . 

The Poisson bracket { , } of two meromorphic functions /, g defined over a sym- 
plectic manifold, is defined by {/, g} :— (V/, JVg); in the Darboux coordinates 
its expression is {/ , g} = Yn=i — §p~§q~- ^ ne P°i sson bracket endows the 

set of first integrals with a structure of Lie algebra. A function F is a first integral 
of (12.11) if and only if {F , H} = (i.e H and F are in involution) . 

A Hamiltonian system with n degrees of freedom, is called meromorphically Li- 
ouville integrable if it possesses n first integrals (including the Hamiltonian) mero- 
morphic over U which are functionally independent and in pairwise involution. 
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2.2. Variational equations. Among the various approaches to the study of 
meromorphic integrability of complex Hamiltonian systems, we choose a Ziglin- 
Morales-Ramis type of approach. Concretely, our starting points are the Morales- 
Ramis MR99] Theorem and its generalization, the Morales-Ramis-Simo Theorem 
[MRRS07] . These two results give necessary conditions for the meromorphic in- 
tegrability of Hamiltonian systems. We need to introduce here the notion of varia- 
tional equation of order m G N along a non punctual integral curve of <|2.1[) . 

Let <fi(z,t) be the flow defined by the equation (|2.1j) . For zq s T, we let 
(j>o{t) := 4>{zq , t) denote a temporal parametrization of a non punctual integral 
curve r of (|2.ip such that zo — (f>(wo,to). We define (VE™ ) the m th variational 

equation of (|2.1I) along T as the differential equation satisfied by the £j := 9 *g z Z j'^ 
for j < m. For instance, (VE| ) is given by (see ; MR99] and [MRRS07]): 



For m = 1, the equation (VE^ q ) is a linear differential equation 

£i = Ai£i where A\ := d^Xjj = J ■ Hess$ {H) <E Sp(n, k) and k := C(0o(*))- 

Higher order variational equations are not linear in general for m > 2. How- 
ever, taking symmetric products, one can give for every (VE?) an equivalent lin- 
ear differential system (LVEJ? ) called the linearized m th variational equation (see 



Since the (LVE™ o ) are linear differential systems, we can consider them under 
the light of differential Galois theory ( |vdPS03),IMR99| ). We take as base field the 
differential field k := C(4>o) generated by the coefficients of cf>o and their derivatives. 
Let K m be a Picard Vessiot extension of (LVE?) for m > 1. The differential Galois 
group G m := Gal(if m /k) of (LVE™ o ) is the group of all differential automorphisms 
of K m that leave the elements of k fixed. 

As G m is isomorphic to a algebraic linear group over C, it is in particular an 
algebraic manifold and we can define its Lie algebra g m := Ti im G° m , the tangent 

space of G m at Id m (with d m = YliLi ("^l^ 1 ) tne s ^ ze °^ (LVE™ )). The Lie algebra 
Q m is a complex vector space of square matrices of size d m whose Lie bracket is 
given by the commutator of matrices [M , N] = M ■ N — N ■ M. We say that g m is 
abelian if [g m , g m ] = 0. 

Following the notations above, we can finally give the Morales-Ramis-Simo 
theorem: 

Theorem 2.1. flMRRS07j ): If the Hamiltonian system 12. 1\) is meromorphi- 
cally Liouville integrable then the g m are abelian for all m G N* . 

Partial effective versions of this theorem have been proposed. In M RRS07] 
(and already [MR99J), a local criterion is given for the case when the first vari- 
ational equation has Weierstrass functions as coefficients ; in |MS09j . a powerful 
approach using certified numerical computations is proposed. In the case of Hamil- 
tonian systems with a homogeneous potential, yet another approach is given in 
[CDMPlOj . 

Our aim is to propose an alternative (algorithmic) method using a (constructive) 
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notion of reduced form for the variational equation. This strategy should sup- 
ply new criteria of non-integrability as well as some kind of "normal form along 
a solution" . We will now explain this notion of reduced form (which we started 
investigating in |AW| ) and show how to apply it. We will then apply our reduction 
method in detail on the well-known degenerated case of the Henon-Heiles system 
proposed in [MRRS07] . 



3. Reduced Forms 

Let (k , ' ) be a differential field with field of constants C and let Y' = AY 
be a linear differential system with A = (ay) 6 A4 n (k). Let G be the differential 
Galois group of this system and the Lie algebra of G. We sometimes use the 
slight notational abuse g = Lie(Y' — AY). 



Let ai,...,a r denote a basis of the C-vector space spanned by the entries 
cti t j 6 k of A. Then we have 



A:=J2 a i( x ) M i> MieM n (C). 



i=l 

This decomposition appears (slightly differently) in jWN63j . we call it a Wei- 
Norman decomposition of A. Although this decomposition is not unique (it depends 
on the choice of the basis (a^)), the C— vector space generated by the Mi is unique. 

Definition 3.1. With these notations, the Lie algebra generated by Mi, . . . , M r 
and their iterated Lie brackets is called the Lie algebra associated to A, and will be 
denoted as Lie(A). 

Example 3.2. Consider the matrix 




Ai := 2(a: 4 -10a: 2 + l) 



Expanding the fraction 2 ^ x x ^.l°^ 1 -jf 1 ^ gives a Wei-Norman decomposition as 






2/x 











2/ 











12a; 
2 + l) a 









where 



Mi = 



Ai = -Mi 

x 



10 

1 

10 





12a; 



-Mi 



M 2 = 







2 

10 



and Lie(Ai) has dimension 6. 

A celebrated theorem of Kovacic (and/or Kolchin) states that g C Lie(A). This 
motivates the following definition : 



Definition 3.3. We say that A is in reduced form if Lie(A) = g. 
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A gauge transformation is a change of variable Y — PZ with P 6 GL(ro , k). 
Then Z' = RZ where R := P^ 1 (AP - P'). In what follows, we adopt the notation 
P[A] := P~ 1 (AP — P') for the system obtained after the gauge transformation 
Y = PZ. 

The following theorem due to Kovacic (and/or Kolchin) ensures the existence of 
a gauge transformation P £ GL(n , k) such that P[A] 6 g(k) when k is a Ci-field3 

Theorem 3.4 fsee |Kov71l lvdPS03j p.25 Corollary 1.32 ). Let k be a differ- 
ential C\-field. Let A G M. n (k) and assume that the differential Galois Group G 
of the system Y' — AY is connected. Let q be the Lie algebra of G. Let H be 
a connected algebraic group such that its Lie algebra f) satisfies A 6 i)(k). Then 
G C H and there exists P G H{k) such that the equivalent differential equation 
F' = AF, with Y = PF and A = P[A] = p- 1 AP - p- l P', satisfies A G g(k). 

We say that a matrix P G GL n (k) is a reduction matrix if P[A] G fl(k), i.e P[A] 
is in reduced form. We say that a matrix Q G GL„ (k) is a partial reduction matrix 
when Q\A\ G f)(k) with g C f) C Lie(A). The general method used to put A in a 
reduced form consists in performing successive partial reductions until a reduced 
form is reached. 

In our paper |AW1 , we provide a reduction algorithm that computes a reduc- 
tion matrix Pi G Sp(2,k) for 4 x 4 linear differential systems Y 1 = A\Y with 
A\ G sp(2,k) (and also for 2 x 2 systems). The first variational equation of a 
Hamiltonian system with n = 2 degrees of freedom belongs to this class of systems. 
If P\ is a reduction matrix for A\ then Sym m P\ is a reduction matrix for sym m A\ 
because Sym m is a group morphism (see [vdPS03 , chapter 2 or FH91 appendix 
B2). 

In what follows, we will assume that we have reduced the first variational equation, 
that its Lie algebra is abelian (so that the Morales- Ramis theorem gives no obstruc- 
tion to integrability) , and use this to start reducing higher variational systems. 

We will follow the philosophy of Kovacic's theorem 13.41 and look for reduction 
matrices inside exp(Lie(A)) . We remark that, in the context of Lie-Vessiot systems, 
an analog of the above Kolchin-Kovacic reduction theorem is given by Blazquez and 
Morales ( BIMolO , section 5, in particular theorems 5.3 and 5.8) in relation to 
Lie reduction. 

The notion of a reduced form is useful in many contexts, such as: inverse 
problems (where the notion was first studied) , the computation of the transcendence 
degree of Picard Vessiot extensions, fast resolution of linear differential systems 
with an abelian Lie algebra and to implement the Wei-Norman method for solving 
linear differential systems with a solvable Lie algebra (using the Campbell-Hausdorff 
formula) WN63 . Reduced forms are also a new and powerful tool that provides 
(non-)abelianity and integrability obstructions for (variational) (see Theorem 12. ip 



A field k is called a Ci-field (or cohomologically trivial) if any homogeneous polynomial 
P £ k[X\ , . . . , X n ]—d of degree d has a non-trivial zero in k n when n > d, i.e the number of 
variables is bigger than the degree. All differential fields of coefficients considered in this article 
will belong to the Ci class. 
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linear differential equations arising from Hamiltonian mechanics, as we will now 
see. 



4. Reduced Forms for Higher Variational Equations 

4.1. Preliminary results. Let (k, ') be a differential field and let d G N. 
Consider a linear differential system Y' = AY whose matrix A G _A/fd(k) is block 
lower triangular as follows: 
(4.1) 



.4 



A 1 

^3 M 



A 



diag' 



-A^,h where A 



diag 



Ax 
A 2 



and A, 




A 3 



The submatrices satisfy A\ G A^di(k), A 2 G A / fd 2 (k), A 3 G Md 2 xd 1 {k) and their 
dimensions add- up d = d\ + d 2 . 



Let 



M 



diag 



Ax 
A 2 



,A t eM dz {k) 



and 




B x 



,Bi G M d2X(il (k) 



Lemma 4.1. Let M X ,M 2 G -M dmg and N 1: N 2 G X sm6 . TTien M X .M 2 e X dl a 9 , 
iVi.iV 2 = (so that Nf = and expf^) = Id + N x ), and [M u Nx] G M sub . 

The proof is a simple linear algebra exercise. 

Let g := Lie(Y' = AY) be the Lie algebra of the Galois group of Y 1 = AY and let 
f) := Lie(A) denote the Lie algebra associated to A. We write t)di ag := rj D Mdiag 
and \) su b ■= fynMsut,- The lemma shows that they are both Lie subalgebras (with 
hsub abelian) and f) = t) diag © t) sub . Furthermore, [f) dia9 , f) sub ] G t) sub (i.e fj^^f, is an 
ideal in f)). When t)diag is abelian, obstructions to the abelianity of f) only lie in 
the brackets [l)di ag , hsub]- 

4.2. A first partial reduction for higher variational equations. Using 
the algorithm of |AWj , we may assume that the first variational equation has been 
put into a reduced form. We further assume that the first variational equation has 
an abelian Lie algebra (so that there is no obstruction to integrability at that level). 

As stated in section [521 each (VE™ o ) is equivalent to a linear differential system 
(LVE™ o ) whose matrix we denote by A m . The structure of the A m is block lower 
triangular , to wit 



(4.2) 



A — 



sym m (Ax) 





Am-l 



eM dm (k) 



where A\ is the matrix of (LVE^ o ). Assume that A m _i has been put in reduced 
form by a reduction matrix P m _i. Then the matrix Q m G GL(d m , k) defined by 



Q, 



Sym m (P 1 ) 






Pm-l 



puts the diagonal blocks of the matrix A m into a reduced form (i.e the system would 
be in reduced form if there were no B m ) and preserves the block lower triangular 



A REDUCTION METHOD FOR HIGHER ORDER VARIATIONAL EQUATIONS 



7 



structure. Indeed, 
where 



Sym m (P 1 )[sym m A 1 ] 

Pm — 1 [-^m — 1 ] 



B m := P m \B m Sym m (P 1 ). 

Applying the notations of the previous section to 1 := Q ra [A m ], we see that 
Lie(A)diag and Lie(A) SU0 are abelian. Obstructions to integrabihty stem from 
brackets between the diagonal and subdiagonal blocks. To aim at a reduced form, 
we need transformations which "remove" as many subdiagonal terms as possible 
while preserving the (already reduced) diagonal part. Recalling Kovacic's theo- 
rem [3Tj our partial reduction matrices will arise as exponentials from subdiagonal 
elements. 

4.3. Reduction tools for higher variational equations. 

Proposition 4.2. Let A := Q m [A m ] as above be the matrix of the m-th 
variational equation Y = AY after reduction of the diagonal part. Write A — 
Adiag + ^2i=i PiBi with Pi e k, where the Bi form a basis of Lie(A) SU0 (in the 
notations of section \4-l\ )- 

Let [A^ag , B\] = X)i=i b liB-i, 7j G k. Assume that the equation y' = jiy + Pi has 
a solution g\ 6 k. Set P := exp(<?i_Bi) = (Ld + g\B\). Then 

d sub 

P[A] = A diag + J2[Pi+ gi ji\ Bi, 



= 2 



i.e P[A] no-longer has any terms in B 



Proof. Recall that P[A] = p- 1 {AP - P') and let P = Ld + g x Bi. We have 
P' = g'\B\ whence 

d Bub 

AP = (A diag + J2 P l B l ){I + gxBx) = A diag + J2 PA + 9iA diag B x 

i=l i>l 

since BiBj = 0. Therefore we have AP — P' = Adiag + g\Adi ag B\ + (Pi - g[)Bi + 
J2i=2 PiBi which implies 

,/ 



P-\AP-P') = (Id-gxBx) 



Adtag + giAdiagBi + {Pi - g[)Bx + ^ P,Bi 



i=2 



= A diag + gi [A diag , B x ] + (Pi - g'^Bx + ^ p i B l 

i-2 

because B x A d iagBx = Bx[A diag , B x \ + A^BxBx = Bi Et^.] = 0. So, as 
9i = 7i5i +Pi, we obtain 

dsiLb 

P[A] = Adiag + £ [Pi + 9lli] Bi. 
i=2 

a 
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Remark 4.3. If 71 = then we simply have gi — J Pi- In that case, suppose 
that k — C(x) and that /?i = R,[ + L\ where R\ € C(x) and L\ € C(x) has only 
simple poles, then J 0i ^ <C(x). However, if we apply proposition 14.21 with the 
change of variable Y = (I + R\B\)Z a term in B\ will be left that will only contain 
simple poles. 

This proposition gives a nice formula for reduction. However, it is hard to i- 
terate unless Lie(A) has additional properties (solvable, nilpotent, etc) because the 
next iteration may "re-introduce" B\ in the matrix (because of the expression of the 
brackets). This proposition provides a reduction strategy when the map [Adiag, ■] 
admits a triangular representation. 

To achieve this, we specialize to the case when the Lie algebra Qdiag has dimen- 
sion (at most) 1. Then we have Adiag = /3qAq where /3o G k and Ao is a constant 
matrix. The above proposition specializes nicely : 

Example 4.4. If A diag = j3 A with j3 G k, A Q e M„(C) and [A , Si] = A.B1 
for some constant eigenvalue A ^ then the change of variable Y = PZ with 
P := (Id + gBi), with g' = \gf3 + ft, satisfies P[A] = /3 A + J2f>2 Pi^i- 

To implement this (and obtain a general reduction method), we let ^0 : \) SU b fysub, 
B i — y [Ao,B]. This is now an endomorphism of a finite dimensional vector space ; 
up to conjugation, we may assume the basis (Bi) to be the basis in which the matrix 
of \l/o is in Jordan form. We are then in position to apply the proposition iteratively 
(see the example below for details on the process). 

Remark 4.5. Not that Aq needs not be diagonal. The calculations of lemma 
14.21 and subsequent proofs remain valid when Aq is block lower triangular. 

We have currently implemented this in Maple for the case when Adiag is 
monogenous, i.e. its associated Lie algebra has dimension 1. We will show the power 
of this method and of the implementation by giving a new proof of non-integrability 
of the degenerate Henon-Heiles system whose first two variational equations are 
abelian but which is not integrable. 

5. A new proof of the non integrability of a degenerate Henon-Heiles 

system 

In this section we consider the following Henon Heiles Hamiltonian [MR99J, 
[MRRS07] . 

(5-1) H := i(p? +p\) + -{ql + q\) + -q\ + \qiql 

as given in [MR 99 . This Hamiltonian's meromorphic non integrability was proved 
in |MRRS07j . The Hamiltonian field is 

q\ = Pi , <72 = p 2 , pi = -qi(l + qi) - ~<?1 , P2 = -92(1 + qi). 

This degenerate Henon Heiles system was an important test case which motivated 
[MRRS07] . Its non integrability was reproved in [MS09 to showcase the method 
used by the authors. We follow in this tradition by giving yet another proof using 
our systematic method. Our reduction provides a kind of "normal form along <p" in 
addition to a non integrability proof. The readers wishing to reproduce the detail 
of the calculations will find a Maple file at the url 
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\protect\vrule widtriOpt\protect\href {http : //www.unilim. f r/pages_perso/jacques-artlrur .weil/ cha 

It contains the commands needed to carry on the reduction of the (LVE™) for 
i = 1 ... 3. The reduction of (LVE^) may take several minutes to complete. 

5.1. Reduction of (VEj,). On the invariant manifold {q2 — , P2 = 0} we 

consider the non punctual particular solution 

3 1 



1,0,4^^^,0 



2cosh(t/2) 2 ' ' 2cosh(i/2) 3 

and the base field is k = C{(/>) = C(e*/ 2 ). Performing the change of independent 
e*/ 2 , we obtain an equivalent system with coefficients in C(x) given by 



variable x 



A x := 






2{x i -Wx 2 + l) 








2/x 











2/x 











12x 
a +1 ja 









Applying the reduction algorithm from |AWj we obtain the reduction matrix 



Pi ■■= 



6(a:-l)(a:+l)3: 2 
(x 2 + l) a 



6x 2 (x i -4x 2 + l) 

(z 2 +l) 4 




E 4 -4a; 2 + l 
(i- 2 + l) 2 



6(x 2 -l)x 2 
(:r 2 + l) 3 



that yields the reduced form 



6x 



z lu + 15z"-16a:''-144a; 4 + 15a.^ + l 
12k 2 (x 2 + l) a 



z 12 +4:r 1 " + 121.x* + 256:r -249a: -4a; 2 - 1 
12a: 2 (a; 2 + l) 4 



10 

6/5 







5a; 4 + 16z 2 -13 
3(a; 2 + l) 2 



x°-x' l -17x- ! +l 
(o- 2 + l) a 



We see that dime (Lie(Ai t n)) = 1 and since ^ has one single pole, we cannot 
further reduce without extending the base field k. We find, 



01 = span c jL>i : 

which is trivially abelian and therefore doesn't give any obstruction to integrability. 



1 

6/5 







5.2. Reduction of (LVE 2 ). We want now to put the matrix A 2 of (LVE 2 ) 
into a reduced form. First we reduce the diagonal blocks as indicated in section 
4.21 using the partial reduction matrix Q2 := Sym a Pl p so that we obtain a 

partially reduced matrix (its diagonal blocks are reduced whereas its subdiagonal 
block is not): 



Q2[A 2 



sym 2 Ai } fi 

B2 A\R 



with 



Q2\A2\diag 
Q2\A2\sub 



sym A ltR 



B 2 



We compute a Wei-Norman decomposition and we obtain an associated Lie algebra 
Lie(Q2[A 2 ]) of dimension 11 such that: 
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On one hand we obtain Lie(Q 2 [A 2 ])di a g — span c ^-D2,o := 
with coefficient flo '■— jj- 

On the other hand, Lie(Q 2 [^2]) sub = s P an c{^2} where 



sym 2 Di 







Di 



10 



B 2 := {Bi := 





Bi 



i = 1 ... 10} and Q 2 [A 2 ] diag = ^ fc,i B %,i with A e k - 



The matrix of the application 

f2,o : Lie{Q 2 [A 2 })sub — > Lie{Q 2 [A 2 ]) sub , Bj i-> [D 2 ,o , Bj] 
expressed in the base B 2 takes the following form: 



*2,o := 







-6/5 











1 




6/5 






We denote by J$, 2 the matrix of 'i , 2,o expressed in its Jordan basis, given by the 



matrices = 
Jordan form is 





C 2 ,i 



and their coefficients 72.4 with i = 1 ... 10. So the 



J\S/, 



0000 



000000 

To perform reduction we will use the Jordan basis C 2 :— {C 2 ,i , i = 1 . . . 10}. The 
decomposition given by the Jordan basis C 2 is Q2[^2] := -Do + Si=i7i^ with 
7i G k , i = 1 . . . 10. We notice that J* is made of three Jordan blocks 
- two blocks of dimension 4 : 



{£2,4 , C 2 ,3 , C*2,2 , £2,1} 



and 



{C'2,8 , C-2,7 
and one block of dimension 2 



C*2,6 1 C*2,5} 
{C2.10 j ^2,9} 

The hypothesis of the first section of Proposition 14.21 are satisfied. Therefore the 
partial reduction of (^[^l is done in the following way: 

- Choose a Jordan block of dimension d : {C 2 ,i ■ ■ ■ C 2 ^+d-i}- It satisfies 

*2,o(C 2 ,i+ s ) = C 2 ,i+ s -i for s = l...d- 1. Set A 2 := Q 2 [A 2 ] and set 

s:=d-l. 
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- For s from d— 1 to 1, compute the decomposition 72,i+s = R'2 i+ s +^2,i+s 
where i?2,i+s , L2,i+s £ k and L2,i+ S has only simple poles. 

Take the change of variable P2A+S = Id + R2,i+ s C2,i+ s and perform the 
gauge transformation i-^i+sl^]- 

If £2,4+5 = then the Wei-Normal decomposition of P2,i+ S [A2] does not 
contain C 2 ,i+ S so C 2 ,i+ S £ 02- 

Set A 2 := P2,i+ S [A2] and set s := s — 1. Repeat this procedure recursively 
until s = 1. 

- Choose a Jordan block that has not been treated. Repeat until there are 
no more Jordan blocks left untreated. 

In this way, only will be left in the subdiagonal block the C2,i that have co- 
efficients L,2 t i (after the procedure) containing only simple poles. In our case, we 
obtain a reduced matrix for (LVE?): Ai,r '■= jC*o and 



C := 



r « 




















12 

3 









Sk 

18 


10 

3 








5 
3 


2 






- 22 

3 



-2 



As in the case of A 1 r, this matrix A 2t n is in a reduced form because Lie(A2 i n) 
is monogenous and \ only has simple poles. Therefore Lie{A2 t n) = 02 and g 2 is 
once more abelian bringing in no obstruction to integrability. We have then to look 
at (LVE»). 



5.3. Reduction of (LVE?). We denote P2 the reduction matrix of A2. Once 



more we build a partial reduction matrix Q3 := 



Sym- s P 1 






Pi 



that puts the 



diagonal blocks of matrix A3 into a reduced form and we obtain the partially 
reduced matrix := sym 6 Al,R . .In this case we have a Wei-Norman 

L J «3 A 2,H 

decomposition of of dimension 18, and dimc(£ie(<33^43])) = 38. 

We thus have 



- dim c (I/ie(Q3[-43])<ija 9 ) = 1 where 

Lie(Q 3 [A 3 ]) diag = span c { J D 3 , := 



Sym-'D 1 
» 62,0 



} 



and dim c (Lie(Q 3 [A 3 ]) sub ) = 37 and Lie(Q 3 [A 3 )) sub = span c (# 3 ) with 



B 3 = {B 3tl = 





B 3,i 



i= 1 ... 38} 



a base of generators of Lie(Q 3 [A 3 ]) su b. 
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We define ^> 3i o : f)3, s «6 — > f)3,s«6 , B n- [Z?3,o , B]. It is nilpotent and its Jordan 
basis will satisfy the conditions of the first section of Proposition l4.2l In the Jordan 
basis C 3 := {63^ , i = 1 . . . 37}, the Jordan form of J$ 3 is formed by the following 
Jordan blocks: 

(1) three Jordan blocks of dimension 5 corresponding to : {C3.5, . . . , £3,1}, 
{63,11: ■ • ■ , Cs.e}, {63,17, . . . , 63,12} 

(2) three Jordan blocks of dimension 4: {6*3,18, . . . , 6*3,21} , {6*3,22, ■ • • , 63,26} 
and {63,31, . . . ,63,27}, 

(3) and two Jordan blocks of dimension 2: 

{63,34, ■ ■ ■ , 63,32} and {63,37, • ■ • , 63,35}. 
In the basis C 3 , a Wei-Norman decomposition is 



37 



Q 3 [^ 3 ] = A)£>3,0 + > 73,i6 3 , 



8=1 



We proceed blockwise as in the case of the second variational equation. This 
time, possible obstructions to integrability appear when handling the Jordan block 
{63,31 , . . . , 63,27}- By decomposition 73, j = R' 3i + L 3yi (with i = 27 ... 31), we see 
that in particular £3,30 and L3 29 are non zero (and have "new poles", i.e not the 
pole zero of the coefficient of the reduced form of (VE2)) and therefore we suspect 
that 63,29,63,30 (or some linear combination) lie in g 3 . Since neither 63,30 nor 
63,29 commute with D 3 that would suggest that g 3 is not abelian and therefore, 
intuitively, the Hamiltonian (|5.ip would be non integrable. We prove this rigorously 
in the following subsection. 

5.4. Proof of non-integrability. After performing the partial reduction re- 
cursively for all blocks, we obtain the matrix A 3 ,r. It has a Wei-Norman de- 
composition A 3 , r — diM 3 ,i + a 2 Af3,2 where M 3 ,i,M 3 ,2 € M. 3 4,(C), cl\ := K 
a 2 '■= x i +l ■ The matrix M 3t \ is lower block triangular and M 3) 2 € Lie(A 3! u) su } ) . 
We let M 3 , 3 := [Af 8 ,i , M 8 , 2 ], M 3A := [M 3 ,i , M 3 , 3 ], M 3 , 5 := [M 3j i , M 3A ) and check 
that [M 3j i , M 3 j] = otherwise. So Lie(A 3t n) has dimension 5 and is generated by 
the M 3) j. Note that M 3jl e M 3 4,. su b{'£) for i > 2. Again we let 

: Lie(A 3 , R ) — * Lie(I 3> fl) , M 1— > [M 3 ,i , M]. 



By construction, the matrix of $ is 



00000 
00000 
01000 
00100 
00010 



Theorem 5.1. ^3,^ is a reduced form for (LVE^) and q 3 is not abelian so the 
degenerate Henon-Heiles Hamiltonian &5.1]) is not meromorphically integrable. 

Proof. We know that Lie(A 3t n) is non abelian so we just need to prove that 
A 3: n is a reduced form. To achieve this we will construct a Picard Vessiot extension 
K 3 still using our "reduction" philosophy and we prove that it has transcendence 
degree 5: as g3 C Lie(A 3t n) and dim.c(Lie(A 3y n)) — 5 this will show that q 3 = 
Lie(A 3yR ) because dim c (fl 3 ) = dtr(if 3 /k) (see lvdPS03] Chap. 1.). 

We apply proposition 14.21 to A 3t R. Apply the partial reduction Pi = {Id + 
J aiM3,i) = Id + ln(a;)M3,i: Pi[A3,/{] contains no terms in M3,2 and -Pi[A3,a] = 
01-^3,1 + («i / 02) ^3,3; we call I2 = J(ai J 0-2) = Li2(x 2 ) where Li 2 denotes the 
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classical dilogarithm (see e.g [Car02]). Similarly we obtain I 3 and I a as coefficients 
of successive changes of variable. We are left with a system Y' = a\MY, the Picard- 
Vessiot extension is 

K 3 = C(x)(ln(x) , ln(l + x 2 ) , Li 2 (x 2 ) , Li 3 (x 2 ) , Li 4 (x 2 )) 

It is known to specialists that dtr(-ff 3 /k) = 5 (and reproved for convenience 
below) . □ 

5.5. A self-contained proof of dtr(i£" 3 /k) = 5. To remain self-contained we 
propose a differential Galois theory proof of the following classical fact (see |Car02j 
for instance). The proof is simple and beautifully consistent with our approach. 
To simplify the notations, we write the proof in the case of the classical iterated 
dilogarithms Lij(-x) but, of course, it applies mutatis mutandis to our case of 
Lij(x 2 ). 

Lemma 5.2. Let K 3 = C(x)(ln(x) , - ln(l - x) , Li 2 (—x) , Li 3 (— x) , Li^—x)), 
then dtr(if 3 /k) = 5 

PROOF. Let us prove that the functions 

x , \n(x) , — ln(l — x) , Li2{—x) , Li 3 (— x) , Lu(— x) 

are algebraically independent using a differential Galois theory argument. That 
ln(a;) and — ln(l — x) arc transcendent and algebraically independent over C(x) is 
a classical easy fact. We focus in proving the transcendence and algebraic indepen- 
dence of Li2{—x) , Li 3 {—x) and Li^—x). Set the following relations, 

Lio(— x) := , Li\{— x) := — ln(l — x), Li 2 (— 2;) := f — — -dx, 



Li 3 (—x) :— [ — — ^-dx, LiA—x) := f — — -dx 

J x J x 

and therefore K 3 = C(x)(ln(x) , Lio(-x), . . . , Li^—x)) is a differential field (with 
Litt-x) = Li *-A- X ) ). Of course, dtr(if 3 /k) < 5. Let us define 

f ln(x) 2 ln(x) 3 

V := span c < 1 , ln(x) , — - — , — - — , Lii(-x) , Li 2 (-x) , Li 3 (-x) , Li&(-x) 



and consider and element a G Gal(K 3 /k). As <r(ln'(x)) = cr(~) — \ — In' (a;) 
there exists a constant cq G C such that <r(ln(x)) = cq. Similarly, we obtain that 
a(ln(x) 2 /2) = ln(x) 2 /2 + c ln(x) + c 2 and cr(ln(x) 3 /6) = ln(x) 3 /6 + c 2 ln(x)/2 + 
Co ln(x) 2 /2cQ. Since Li[(—x) = xi x +1 G k we have that a(Li[(—x)) — Li'^—x) and 
therefore there exists c\ G C such that <j{Li\{—x)) = Lii(—x) + c\. As Li' 2 (—x) = 
£tl ^ 1} we have that a(Li' 2 (-x)) = a( £tl ^ 1] ) = Ul( x ' x) + f and there exists 
c 2 G C such that ct(Li 2 (— x)) = Li 2 (—x) + Cim(x) + c 2 . We prove similarly the 
existence of c 3 , C4 G C such that 

ln(x) 2 

a(Li 3 (-x)) = Li 3 (~x) + c 1 — - h c 2 ln(x) + c 3 

t-r- 1 \ \ T ■ , \ hi(x) 3 ln(x) 2 , 
<7{Li A {-x)) = Li4,(-x) + ci — hc 2 — h c 3 m(x) + c 4 . 
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We see that V is stable under the action of Gal{K$/\s) and hence is the solution 
space of a differential operator L G k[^] of order 8. Therefore, in this basis the 
matrix of the action of a on V: 



M a := 



c 

1 










=o/2 
CO 

1 









3/6 



ci 




1 







C2 
Cl 




1 






C3 
C2 
Cl 



1 
H 



C4 
C3 



As \vl(x) and ln(l — x) are transcendent (and algebraically independent) we know 
that cq and c\ span C. It follows that 33 contains at least 



mo := 



and toi 



Since too and Toi do not commute, we know that the Lie algebra generated by the 
iterated Lie brackets has dimension at least 3. Iterating the brackets of Too and toi 
we obtain a subalgebra of q% of dimension 5. Therefore we have dtr(ii3/k) > 5 and 
since we know that dtr(X3/k) < 5 we obtain the equality and the result follows. 

□ 



Remark 5.3. Horozov and Stoyanova [HS07 make use of the properties of the 
dilogarithm in order to prove the non-integrability of some subfamilies of Painleve 
VI equations: namely, they prove the non-abelianity of Q2, the Lie algebra of its 
second variational equation. 



6. Conclusion 

The reduction method proposed here is systematic (and we have implemented 
it in Maple). Although it is currently limited to the case when Lie((VEj,)) is one- 
dimensional, extensions to higher dimensional cases along the same guidelines are 
in progress and will appear in subsequent work. In work in progress with S. Simon, 
we will show another use of reduced forms, namely the expression of taylor expan- 
sions of first integrals along (j) are then greatly simplified. 

We conjecture that our method is not only a partial reduction procedure but a 
complete reduction algorithm : assuming that (LVE^ ) is reduced (with an abelian 

Lie algebra), we believe that the output A m+ i.R of our reduction procedure of sec- 
tions 4 and 5 will always be a reduced form. In the context of complex Hamiltonian 
systems, this would mean that our method would lead to an effective version of the 
Morales-Ramis-Simo theorem. 
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